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Abstract 

The response of a field-stabilized two-dimensional smectic to shear stress 
is discussed. Below a critical temperature the smectic film exhibits elastic 
response to an infinitesimal shear stress normal to the layering. At finite 
stresses free dislocations nucleate and relax the applied stress. The coupling 
of the dislocation current to the stress results in non-newtonian viscous flow. 
The flow profile in a channel geometry is shown to change qualitatively from a 
power-law dependence to a Poiseuille-like profile opon increasing the pressure 
head. 

PACS numbers: 68.10.Et, 47.50.+d, 83.70.Jr 

I. INTRODUCTION 

Thin layers of elongated molecules exhibit a variety of different liquid crystalline phases 
differing in the degree of positional and orientational order from isotropic films. Complex 
rheological properties are common, in particular, when external orienting fields affect the 
alignment of the molecules. The coupling of the velocity field to the liquid-crystalline or- 
der gives rise to anisotropic viscosities. Furthermore, nonlinearities introduce shear-rate 
dependent effective viscosities, which includes phenomena like shear thinning, etc. 

Recent experiments on Langmuir monolayers []!]-§[ on an air-water interface demonstrate 
the flow-induced microstructural deformations and orientation within a liquid- crystalline 
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film. In particular, the non-linear shear response gives rise to unconventional flow profiles 
in channel geometry. Since dislocations often play a prominent role in two-dimensional 
systems, the motion of these point defects is presumably responsible for many of the observed 
phenomena. 

Here we focus on a smectic film floating on a liquid substrate. Smectic films are two- 
dimensional layered systems characterized by an elastic free energy that descibes long- 
wavelength distortions of the layers in terms of a displacement field u(r) [|J. Up to quadratic 
order, we have 

F= 1 -Jd*r[B(d z u)* + B\*(dlu) 2 }, (1) 

where the first term deals with layer compression, while the second term is related to splay 
distortions. The thermal fluctuations destroy even quasi-long-range translational order at 
any nonzero temperature. One easily finds that phonons give rise to exponential decay in 
the order parameter ip(r) = exp(iq u(r)) , where go — 27r/<i is the wavevector associated with 
the layer spacing d []5| . In order to moderate the effect of thermal fluctuations one can apply 
an in-plane orienting external field. The corresponding free energy then has to be modified 
to § 

F=\j d 2 r[B{d z uf + E(d x u) 2 } . (2) 

Here E measures the coupling to the external field, e.g. an electric or magnetic field. 
The nematic director orientation is aligned with the displacement field according to 9 = 
—const. d x u on large scales. This property does not necessarily hold near the boundaries. 
As discussed in the Appendix, there appears a characteristic length scale on which the 
director adjusts to the displacement. 

One can eliminate the anisotropy in Eq. (|2|) by a simple volume-preserving rescaling, and 
the static properties are thus determined by the universality class of the AF-model. There is 
a critical temperature T c oc V BE separating a low-temperature smectic phase characterized 
by bound dislocations and power-law correlations in (tfj(r)*tfj(0)) and a high-temperature 
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nematic phase with exponential decaying order parameter correlations and free dislocations. 
Varying the external field allows control of the anisotropy. In particular, since the ratio E/B 
is preserved under renormalization, one can adjust experimentally the critical temperature 



Here we discuss the dynamical properties of the stabilized smectic below T c . The long- 
wavelength renormalized stiffness constants Er and Br now attain finite non-zero values and 
the film reacts as an elastic medium to an infinitesimal applied shear stress perpendicaular 
to the layers. However, a finite stress can trigger the proliferation of free dislocations which 
results in shear flow. 



The long-wavelength and low-frequency dynamics of a smectic film is governed by hy- 
drodynamic equations for the broken symmetry variable, viz. the layer displacement u, 
and the conserved quantities Q. For good thermal conductivity on a water substrate no 
temperature gradients can build up and energy conservation can be ignored. Furthermore 
we specialize to incompressible smectic films, setting the mass density p = const. Mass 
conservation then implies the constraint d x g x + d z g z = for the momentum density g(x, z). 
Similarly the trace of the viscous stress tensor also drops out of this long-wavelength low- 
frequency description. The viscous stress tensor then has only two independent components 
which have to be related to the two of the strain rate tensor. From symmetry there appear 
two indpendent viscosities and we write the constitutive equation for the viscous stress as 
a 'xz = u i.9 z g x + d x g z ),a' xx — {y — v')d x g x ,a' zz — (v — u')d z g z . The linearized equations of 
motion then read 



T 
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II. HYDRODYNAMICS 




(3a) 



dt9. 



X 



d x p + vd\ 



9x + v'd x d z g. 



(3b) 



d t g z = -d z p + {Bd 2 z + Ed 2 x )u + vd 2 x g z + u'd z d x g,. 



X ) 



(3c) 
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where X p denotes the permeation constant and p the pressure. 

If a constant external shear stress is applied, the medium reacts by building up a 
stationary viscous momentum flow o e x x ^ = vd z g x . Since the layers are not distorted no elastic 
stresses can build up to balance the external stress. However, the medium does sustain an 
infinitesimal external stress cr||* by straining layers according to o e * x = Ed x u. 

The preceding paragraphs correctly describe the linear response of the stabilized smectic 
film. The nonlinear response, however, is qualitatively different. Finite applied stresses 
create free dislocations in the displacement field which can move and relax the stress. For 
two-dimensional smectics dislocations have Burgers' vector ±d along the z-axis (see Fig. 
1). Hence, we can think of these dislocations as scalar defects with charge m v = ±1. In 
the presence of plastic flow due to the dislocation motion the film now exhibits viscous 
behavior. A convenient way to incorporate the effects of dislocations is by switching from 
a description in terms of the displacement field u to the strains s x ,s z . This avoids the 
introduction of branch cuts in the displacement field u(x,z). Locally the strains are given 
by s x = d x u, s z = d z u. However, since the displacement field is no longer a single-valued 
function, the line integral § T {s x dx + s z dz) does not vanish for closed loops T, but rather 
counts the number of enclosed dislocations in units of the layer spacing d. Consequently the 
curl of the strain is given by d x s z — d z s x = m(r)d, where m(r) = J2v=±i rrivn v {r) is the total 
"charge" dislocation density. The charges m u = ±1 characterize the point dislocations and 
n v (r) = $(f—r iv ) is the number density of dislocations of charge m v . Since dislocations 
can be created only in pairs, the overall dislocation density is conserved, 



Here J denotes the two-dimensional current of dislocation charge. In terms of the strains, Eq. 
(ft), is merely the integrability condition which allows us to set d t s z + J x d = d z E, d t s x — J z d = 
<9 X 5. Choosing H = g z j p+X p (Bd z s z + Ed x s x ) ensures compatability with the dislocation-free 



d t m + div J = . 



(4) 





(5a) 
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d t s z = d z g z /p + \ p d z (Bd z s z + Ed x s x ) - J x d , 



(5b) 



d t g x = -d x p + vd 2 z g x + v'd x d z g. 



z , 



(5c) 



d t g z = -d z p + {Bd z s z + Ed x s x ) + vd 2 x g z + v'd z d x g. 



(5d) 



For the dislocation current we adopt a Fokker-Planck description. There is a diffusive 
current due to gradients in the dislocation density. Symmetry dictates that the principal axis 
of the diffusion tensor are aligned with the (x, z)-coordinate system of the layers. Further- 
more stresses induce motions that result in a separation of dislocations of opposite charge. 
Following Ref. p we write 



The first terms on the right-hand side are known as the Peach-Koehler forces and are the 
analog of the Magnus force in a superconductor, see Fig 1. The diffusion constants D x = 
r x ^T and D z = T z ksT refer to dislocation climb and glide. Since climb, i.e. motion 
perpendicular to the Burgers vector, involves long-range mass transport one expects T x <C 
r z . The free dislocation density n(r) = ^ v Tn 2 u n v {r) vanishes in the equilibrium smectic 
phase since all dislocations are bound. Thus the linear dynamics given by Eqs. (|5|) and (^|) 
merely adds an anisotropic diffusive mode for dislocation motion to Eqs. (Q). 

A non-zero free dislocation density n(r) leads to viscous response to an external shear o e zx l 
as shown in Ref. . Shearing perpendicularly to the layers of the smectic film to produce a 
nonzero s x ~ d x u liberates free dislocations which then can relax the shear. Thus one has a 
most prominent form of the phenomenon known as shear thinning, i.e. a reduction from an 
infinite linear viscosity and a non-zero linear elastic modulus to a finite nonlinear viscosity 
and zero nonlinear elastic shear modulus. This is the subject of the next two sections where 
we adopt methods applied to superfluid helium films in Ref. to smectic hydrodynamics. 



Jx = T x Bds z n — T x kBTd x m , 



(6a) 



T z Eds x n — r z kBTd z m . 



(6b) 
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III. NUCLEATION RATE OF FREE DISLOCATIONS 



The free energy of Eq. (EI) contains contributions from the smooth part of the strain 
and from the interaction energy of dislocations. The strains can be decomposed as s x = 
d x (j) + s x ng , s z = d z (p + s s z mg . Then the smooth part and the dislocation interaction decouple 
provided one imposes the stress equlibrium condition 

Ed x s x ing + Bd z s s z m9 = . (7) 

This relation is in fact the integrability condition which allows the introduction of the Airy 
stress function Es x ng = —d z x, Bs s z ng = d x x- Since the curl of the strain yields the disloca- 
tion density, the Airy stress function is given in terms of the solution of 

1 1 

-g d lx + -£plx = m{r)d . (8) 

After a partial integration the free energy, Eq. (0), can be written 

T = F - - J d 2 rx{r)m(r)d , (9) 

where the smooth fluctuations of the displacement field are encoded in 

T G = \ J d 2 r[B(d z 4>) 2 + E(d x 4>) 2 ] . (10) 

The dislocation contribution is the smectic analog of the interaction energy of charges in 
electrostatics. Here, m(f) corresponds to the electric charge density and x{^) to the electro- 
static potential. As mentioned in Sec. H the electrostatic analogy becomes complete after 
the volume prserving transformation x = a£,z = (/a, with a = (E/B) 1 ^. Equation @ 
then reads 

dfx + d 2 x = mdBa 2 . (11) 

The bare interaction energy of a pair of dislocations of opposite charge ±1 is then easily 
calculated to be 
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U (g) = HQ/ao) + 2E c (a ) , (12) 

where g = V£ 2 + C 2 is the distance of the dislocations in scaled units and ao ~ d is a short 
distance cutoff that signals the breakdown of continuum elasticity theory. The core energy 
2E c (clq) plays the role of a chemical potential, i.e. the energy to create a defect pair at 
distance a relative to a dislocation free system. Since there are many pairs the effective 
interaction of a particular pair is a many-body problem. The pair carries a polarization cloud 
of bound pairs which screen the bare interaction. Following Kosterlitz and Thouless P this 
problem can be dealt with by the introduction of a scale-dependent dielectric constant e(g). 
The effective interaction then reads 



TT , , Ba 2 d 2 rs dp' . > , . 



The probability T(^) per unit area to find a pair at distance g is then determined by the 
effective interaction T(g) = a^ 4 exp(— ?7(^)//c B T), or equivalently 

(TT = Ba 2 d 2 

2nkBTe(g)g 

In terms of a scale dependent stiffness K(g) = Ba 2 d 2 /[4:7i 2 kBTe(g)} and the dislocation 
fugacity y, where y(g) 2 = g 4 T(g), this is the first of the celebrated Kosterlitz recursion 
relations 

^ = [2-^b. (15) 

The renormalization of the stiffness is due to polarization of the dislocation pairs and is 
governed by the second recursion relation 

dK~ l 



47rV . (16) 



ding 

The flow equations for the dielectric constant and the fugacity then reveal the existence of a 
low-temperature phase with a finite long-wavelength dielectric constant. Since e(g) depends 
intrinsically only on the flow parameter / = ln(^/ao) the first term in the effective potential 
is basically logarithmic and binds the pair only weakly. 



In the presence of a strain s x the pair is subjected to the Peach-Koehler force Es x d that 
separates the pair in the ^-direction in addition to the attraction force corresponding to Eq. 
(|T3|). For non-zero s x the total potential exhibits a saddle point on the z-axis and the pair 
can break up by escaping over the barrier as illustrated in Fig 2. This creates a pair of free 
dislocations that contributes to the relaxation of the applied shear. 

The motion of bound pairs is again given in terms of a Fokker-Planck equation . Since 
climb motion is much slower than glide T x <C T z the associated currents are predominately 
unidirectional. With this simplification one derives 

d t T = d z J z , 

Jz = -2T z T[d z U + Es x d] - 2k B TT z d z T . (17) 

Since we are dealing with the relative motion of a pair, the diffusion constant is 2D Z , where 
D z = T z ksT refers to diffusion of a single dislocation. The escape rate of this effectively one- 
dimensional problem can be obtained by standard methods. The saddle point Zq = (o/ a of 
the total potential determined by the implicit relation d+27ie(\( \)(oS x a = 0. For definiteness 
we discuss the case s x < 0, i.e. the saddle point lies on the positive z-axis. The Fokker- 
Planck equation QT7| ) implies that J z is independent of the z-coordinate. The solution can 
be obtained in terms of if) — [U + Es x zd]/kBT by writing 

J z J e^dz' = -2k B TT z e^r\™ . (18) 

For z = +oo we expect that the dislocation probability density T vanishes rapidly, whereas 
for \z\ <C Zq it assumes its equlibrium value. The right-hand side of Eq. ([18]) then simplifies 
to 2kBTT z a ~ 4 . The left-hand side is dominated by the saddle point and the integral can 
be extended to — oo. Upon expanding near the saddle point and using the saddle-point 
equation, the Boltzmann weight can be approximated by 

M -27rX(Co) r 2 2/ _ \2 

'^W-'^^^-^^V^ 1 ' (19) 

Here we neglected terms involving the derivative of the dielectric constant, since these are 
small on large scales according to Eq. fll6D- The dislocation current is then found as 
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Jz = Ry K((o) /(a( Q ) ex-p[—TcK(( )x 2 /(a 2 ^)], where the production rate R = f^J'zdx 
reads 

R = 2k B TT z a 2 y{C,fQ i e 2 ^ . (20) 

For small strain s x and not too close to the critical temperature, one finds £o ^> where £_ 
is the correlation length implied by the recursion relations, Eqs. (|H||T6l) . Then K(( ) can be 
safely replaced by its large distance limit K(g = oo) = 2/tt[1 +x(T)/4]. Here x(T) measures 
the distance to the critical point, and for temperatures close to T c one has x(T) ~ (T c — T) 1 / 2 . 
Equation (15) then implies in this regime y(g) ~ g- x ( T )/ 2 and correspondingly the production 



rate exhibits a power-law dependence on the strain 

R~ \s x \ 4+x{T) - (21) 

The dynamics of the free dislocation density is governed by the rate equation 

d t n = R - rn 2 . (22) 

The recombination process of two free dislocations of opposite charge is encoded in the rate 
constant r. Since local equlibrium is established much faster than dynamics of the broken 
and conserved variables, one can assume that the free dislocation density follows the rate 
adiabatically, i.e, n ~ R}l 2 ~ \s x \ 2+x ^ T ^ 2 . 

A stationary, uniform external stress a e z xt is balanced by a uniform strain Es x , Eq. 
(jSHp. This results in a nonlinear dislocation flow J z ~ sj, 7 = 3 + x(T)/2, which gives 
rise to shear strain rates d x g z = —J z pd. Upon collecting terms one derives the nonlinear 
constitutive equation 

d x g z ^{a e z f)\ (23a) 
7 = 3 + x(T)/2. (23b) 
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IV. CHANNEL FLOW 



In this section we discuss the application of these ideas to flow of the stabilized smectic 
due to a pressure gradient along a channel. The film layers are oriented perpendicularly to 
the channel flow, i.e. the pressure gradient is along the z-direction. Furthermore we assume 
that the channel is much wider than the correlation length £_ so that the hydrodynamic 
description is valid. 

In the steady state all z-derivatives vanish due to translational invariance along the 
channel, except for the pressure head n' = —d z p = const. The equations of motion, Eqs. 
(|), then reduce to 

= d x g z /p + \ p Ed 2 x s x + J z d , (24a) 

= tt' + Ed x s x + vd 2 x g z . (24b) 

Furthermore we impose the no-slip boundary condition g z (±a) = and vanishing permeation 
current d 2 s x {±a) = 0, where ±a are the walls of the channel. Symmetry then implies 
d x g z (0) = d 2 s x (0) = s x (0) = 0. The derivative d x g z can now be eliminated from Eq. (|24"D 
yielding 

- n'x/E = s x - 5 2 d 2 x s x - puJ z d/E , (25) 

where we introduced the permeation length 5 = (pv\ p ) l l 2 . Since the dislocation current 
is nonlinear in the strain, the flow profile depends on the pressure head. The combination 
pvJ z djE can be written as — s x |As x | 7_1 , where A is a dimensionless constant. Rescaling of 
s x reveals another characteristic length scale L n = E / (Att') , apart from the channel width 
and the permeation length. There are several cases that can be discussed analytically. 

For L n ^> a the pressure head is mostly balanced by elastic deformations of the smec- 
tic and the strain is given approximately by s x = —ir'x/E. The dislocation current is 
negligible compared to the elastic contributions, however, it gives rise to the fluid flow 
d x g z = Es x \As x \' y ~ 1 jv. The flow profile then has the ineresting non- Newtonian form 
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where 7 is a continuously variable exponent in the range < 7 = 3 + x(T)/2 < 3. In 
particular, the velocity is non-linear in the pressure head contrary to conventional Poiseuille 
flow. 

In the opposite case L n <C a, we assume that the permeation current can be neglected, 
except for a boundary layer of order 5. Then for \x\ <C L n , one finds again s x = —nx/E 
whereas for \x\ ^> L n one observes Poiseuille behavior d x g z = —ir'x/v. The solution is 
self-consistent provided 5 <C a(a/ L^ 7-1 ^^ . 

V. CONCLUSION 

We derived the hydrodynamic equations of motion for a smctic film stabilized by an 
external field. Below the Kosterlitz-Thouless transition dislocations are bound in pairs. 
However, a finite shear stress applied perpendicularly to the smectic layers nucleates dislo- 
cation pairs. The motion of these free dislocations results in a non-linear viscous response. 
The smectic film exhibits shear-thinning giving rise to unusual flow profiles in channel flow. 
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APPENDIX A: ALIGNMENT OF THE NEMATIC DIRECTOR WITH THE 

DISPLACEMENT FIELD 



The free energy Eq. (0) is the long-wavelength description of the smectic film. The 
nematic director 9 is locked to the displacement via 6 = —const. d x u. This can be seen 
by taking a more microscopic approach by a free energy that includes the nematic director 
degrees of freedom explicitly, f||j 

F=^J d 2 r[B(d z u) 2 + D(d x u + 6) 2 + K (V6) 2 + E6 2 } . (Al) 

Here D measures the free energy cost for nonaligned director to displacement, Kq is the 
Franck constant and E the coupling of the director to the external field. For an infinite 
system the contribution the director field can be integrated out leading to an effective free 
energy in terms of the Fourier transformed displacement field 

'=U™ + ^- E + £> 4 + K# Mi >r- (A2) 

Expanding to lowest nontrivial order in powers of the wavevector q leads to the free energy 
Eq. (D with E = DE/(D + E). 

Near the boundary we use the more microscopic free energy Eq. ( |A1| ) and derive the 
corresponding Euler equations 

Dd x u - KV 2 9 + (D + E)6 = . (A3a) 
Bd 2 u + Dd 2 x u - Dd x 6 = , (A3b) 



The first equation reveals the characteristic length scale Ao = y KJ (D + E). For variations 
of the displacement field with wavevectors q Aq 1 the director field is given by 6 = 
—d x iiD/(D+E). Near the boundary where independent boundary conditions on the director 
and the displacement field can be imposed is a layer of order Ao where the previous relation 
does not hold. The second equations shows that the parameter entering the Euler equations 
corresponding to the effective free energy is E = DE/(D + E). 
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FIGURES 

FIG. 1. Layering of the elongated molecules of a two-dimensional smectic. The dislocation on 
the right half has positive "charge" , whereas the dislocation on the left possesses negative "charge" . 
The strain s x exerts a Peach-Koehler force on the free dislocations indicated by the arrows. The 
induced dislocation current then relaxes the strain. 

FIG. 2. Sketch of the effective potential for a dislocation pair in the presence of shear 
ksT^p = U + Es x zd. The heavy line marks the boundary of bound and free dislocation pairs. 
Breaking of a bound pair occurs via escape over the saddle point. Since climb motion is negligible 
compared to glide the x-coordinate is effectivley frozen. 
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